ON 



JLAB-THY-97-08 

WM-97-106 

CFNUL-97-02 



Effects of Relativistic Dynamics in pp — ► ppir near Threshold 



J. Adam, Jr. li2 , Alfred Stadler 3 , M. T. Pena 3 ' 4 , Franz Gross 1,5 

1 Jefferson Lab, 12000 Jefferson Avenue, Newport News, VA 23606 

2 Institute of Nuclear Physics, Rez n. Prague, CZ-25068, Czech Republic 

3 Centro de Fisica Nuclear da Universidade de Lisboa, 1699 Lisboa Codex, Portugal 

4 Departamento de Fisica, Instituto Superior Tecnico, 1096 Lisboa Codex, Portugal 

5 Department of Physics, College of William and Mary, Williamsburg, Virginia 23185 

(February 9, 2008) 

<N . Abstract 

> . 
I/"") \ The cross-section for threshold 7r° production in proton-proton collisions 

is evaluated in the framework of the covariant spectator description. The 
{Sj ■ negative energy intermediate states are included non-perturbatively and seen 

to yield a considerably smaller contribution, when compared to perturbative 
q^ , treatments. A family of OBE- models with different off-shell scalar coupling 

is considered. 
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I. INTRODUCTION 

The reaction pp — ► ppn° near threshold is very sensitive to the description of the 
iVTV-interaction and to 7r° production mechanisms. Due to the underlying chiral symmetry 
both the single-nucleon term and the pion-rescattering contribution are suppressed near 
threshold [0-0]. The transition amplitude then results from a delicate interference between 
these and various additional contributions of shorter range [@,f|-10]. The relevance and 



relative importance of various reaction mechanisms, in particular the short range ones, has 
not yet been firmly established. A treatment of these mechanisms consistent with the NN 
interactions employed in the description of the final and initial state is needed to clarify the 
large model dependence. Moreover, the magnitude of the transferred hadronic momentum, 
which is typically p/m ~ J fi/m ~ 0.4, calls for a relativistic approach. 

In this letter we re-examine the single-nucleon term using the covariant spectator de- 



scription [12- 1 I] for both the production mechanism and for the initial and final state pp 



interaction. The pion emission from a single proton is described by the Feynman diagrams 



of Figure 1. The OBE spectator models [13| allow a direct covariant evaluation of diagrams 
(a)-(c). Dynamically, these are the only model-independent contributions. There might be 
additional reaction mechanisms, in particular involving pion rescattering, such as diagram 
(d). They have not yet been included in our calculations. The intermediate state nucleons 
in diagrams (a)-(c) can propagate with positive and negative energy. We will split the con- 
tribution of these diagrams into amplitudes a ± ,b ± and c ' with the superscript indicating 
the sign of the energy of the off-shell nucleon after and/or before the pion emission. 

Our approach differs crucially from earlier ones in two aspects. Firstly, it includes non- 
perturbatively diagrams with intermediate negative-energy nucleons. In perturbative ap- 
proaches these contributions are simulated by the inclusion of effective meson-exchange 
operators acting in two-nucleon space (their contributions are often called "Z-diagrams"). 
More specifically, one takes the amplitudes a~ and b~ (with the nucleons in the initial and 
final state on-shell) to the lowest order in v/c ~ p/m, replacing the NN interaction by 
a single one-meson exchange, usually taking into account only the most important meson 
exchanges. Written in the space of two-component Pauli spinors these amplitudes are iden- 
tified with effective pion-production meson-exchange operators and their expectation values 
are then calculated with conventional nuclear wave functions. These transition amplitudes 
simulate those contributions of the diagrams a~ ,b~ , c H and c h for which the transition to 
negative-energy nucleons occurs only once. 



In previous perturbative calculations [f2|-|T(H] the most important contributions to the Z- 
diagram operator were found to be the exchanges of a, u, and n mesons. It has been shown 
that the pion Z-graph yields a much smaller contribution than the a and uo diagrams. It 
is often considered together with other model dependent off-shell pion-rescattering processes 
PJ9|JT0H . In t>/c-expansion calculations, the a and uj Z-graphs were found to dominate the 
short range reaction mechanisms and provide contributions essential for the explanation of 
the data. In contrast, in the present covariant calculations we obtained significantly smaller 
contributions from terms with negative-energy nucleons. 



Secondly, the NN- models considered recently [i I] include non- minimal dynamical effects 



produced by off-shell scalar-nucleon coupling, scaled by a parameter v. 



T s (p > ,p)=g s l [ l + ^- [/+^-2m]| , (1) 

with v a = —0.75 v, and v$ = 2.6 is. Both iVW-observables and the triton binding energy 
appear to be rather sensitive to this off-shell extension, hence it is very interesting to find out 
if it also works favorably in the considered reaction. Taking only diagrams (a) and (b) for o 
and uj exchange in Born approximation, we, indeed, have found a rather strong dependence 
on the parameter is, mainly due to the cancellation between a and oj terms. However, the 
variation is much weaker when the full interaction is included. 



II. THEORETICAL DESCRIPTION 

Near threshold the reaction is dominated by transition between initial 3 Pq and final 1 S'o 
pp partial waves; only these are considered in our calculations. The total cross section is 



given |15j in terms of the covariant transition amplitude Ai: 

p(T lab ) , (2) 

a 2 M* x p, x ,\^\ 2 , (3) 

where m is the nucleon mass; p,E = y/p 2 + m 2 are the relative momentum and energy in 
the initial state; p(Ti ab ) is the phase space density, T iah is the kinetic energy of the incoming 
proton. The amplitude M. and the phase space density p have dimensions MeV -3 , and MeV 4 , 
respectively. The helicity transition amplitude -Mvai,a'a 2 is calculated from the Feynman 
rules as defined in flBjj. For the partial waves considered, there is only one independent 



a[pb] = 


_ 10 4 (M 2 m 2 TXW 
(2vr) 5 4pE ' ' 
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helicity amplitude, M. J + + ++ = i exp(i 63(E)) Ai, with 63(E) being the 3 Pq phase shift 



The iViV-interactions |]T3[ used in our calculations were fitted to np-data and the effect 



of the Coulomb interaction is not included. It would suppress the cross section at threshold 
by about 40%. It has been shown [|5], that using a constant amplitude in (|2|) results in a 
roughly quadratic dependence of a on T^, clearly at variance with the experimental data 
(Fig. §. 

The most pronounced energy dependence of the total transition amplitude comes from 
the energy dependence of the final state 1 Sq interaction. We deal with it in an approximate 
way similar to that employed in |T^]: we assume that the energy dependence of the final 



state scattering matrix is well approximated by the energy dependence of the corresponding 
on-shell NN scattering matrix M(p) = M(p,p), 

Expressing M(p') in terms of the 1 5'o phase-shift 6± and M (0) in terms of the singlet scat- 
tering length a s , one obtains 

M(p',p) « -ex P ( i 6 1 (E'))^ Si ^ 6l{E ' )) M(0,p), (5) 



with p', E' being the relative momentum and energy of the final pp state. Accordingly, the 
square of the transition amplitude contains the energy-dependent factor 

ct >\ ^ 2 sin 2 (<5i(E')) 

^ ) = ( ^ } {a s p>y ■ (6) 

We introduce it into (|2|) by replacing 

,_ , ,. _ , f m , m . dq 

p{Tiab) — ► p{$, } T lab ) = J —dp 1 —dp 2 — 

f(p') 5 {i \pi+p 2 -p[-p' 1 -q). (7) 

The difference between p(l, T/ a b) and p(£, Ti a b) is illustrated in Fig. ||. The approximation 
introduced in ([J) implies for the Heitler i?-matrix R(p',p)/R(p') « R(0,p)/R(0). Figure 
3 demonstrates the quality of this assumption about the energy dependence of the off- 
shell .R-matrices. The curves corresponding to different energies would coincide in case the 
factorization describing the energy dependence of the half-off-shell NN amplitudes R ++ (p',p) 
were exact. 

We note that we did not employ the effective range approximation for the final state 
scattering amplitude as it was used in Ref . |3| ; it was found in Ref . || to artificially increase 
the cross section. Since the approximation (|4j) allows for a very good description of the energy 
dependence of the data, it can be used to extract the value of the threshold amplitude M.. 
A fit to the data of Ref. [|IBJ yields a value 



\M exv \ = 30.5 x KT 7 MeV- J , (8) 

which is to be compared with our calculations summarized in Table |. A similar fit employing 
the unmodified phase space density with an unrealistic energy dependence of the cross section 
leads to a much smaller value |-M exp | = 5.0 x 10 -7 MeV~ 3 . 

III. RESULTS AND CONCLUSION 

Table | shows in detail the contributions of diagrams (a)-(c) for several models that differ 
in their off-shell parameter v. The best fit to the NN data was obtained ]n| by models 
with v in the range 1.5 — 2.0. From the models listed in the table, the one with v = 1.6 



yields the triton binding energy closest to its experimental value [14]. 

The amplitudes are, in general, complex numbers, but we can subtract the two contri- 
butions from diagram (c) which correspond to having both intermediate nucleons on-shell 
either before or after pion is emitted, and add them to the amplitudes of diagrams (a) and 
(b), respectively. Then we arrive at new partial amplitudes that can be written as real 
numbers (which are given in Table |) times a common phase factor. However, one has to be 
cautious when comparing to other calculations that present results for the original diagrams. 
At threshold, the final state phase shift is zero and in terms of the real numbers A, B, and 
C the actual partial amplitudes are given as 



a ± = A ± / cos 83 

b~ = exp(iS 3 ) B~ 
c ±,+ = exp(i5 3 ) (C ±,+ + iA ± tanS 3 ) 
c ±,_ = exp(i<5 3 ) C ± '" . (9) 

The upper part of the table presents results of our covariant calculations. Diagram (c) 
dominates at threshold. Our results show a moderate dependence on the parameter v and 
we observe that the amplitudes for the models with nonzero v are somewhat enhanced. The 
largest value, at v = 1, is about 20 % larger than the smallest at v = 0. This translates 
into a more than 40 % variation in the total cross section. Although in itself it represents 
a sizeable effect, it is at this point not clear how much of it will survive in more complete 
calculations that would include other, model-dependent reaction mechanisms (such as that 
of diagram (d)). 

Our total result is 2-3 times smaller than the fitted amplitude (|||) for all models. The 
relativistic NN potentials |Tj| contain form factors for off-shell nucleons (sometimes called 
"sideways" form factors). For consistency we include them also in the pion emission vertex. 
They decrease the results by up to 10 % (see Table |). For reference and for comparison with 
perturbative results the total results without these form factors are also given (line " no N 
ffs "). The nucleon form factors are not included in our perturbative transition operators. 
It is difficult to take them into account consistently in v/c- expansion approaches, and the 
usual NN potentials and pion emission operators never consider them. 

The last three lines of Table | list the usual perturbative results obtained with the same 
nucleon wave functions and with the meson parameters and meson-nucleon form factors 
consistent with the A^ A" interaction employed. The first line gives the perturbative IA, to be 
compared with the subtotal M. + = A + + C ++ of covariant results. The perturbative result 
is somewhat larger, the difference between these numbers is due to the w/c-decomposition 
in the one-nucleon operator at threshold (in the covariant description it contains a factor 
1/E compared to 1/m in the perturbative one) and to the presence of the off-shell nucleon 
form factors in the relativistic calculations. 

The perturbative contribution of the Z-diagrams given in the second last line is to be 
compared to the sum of all terms with negative energy nucleon(s) in intermediate states: 
M.~ = A~ + B~ + C + ~ + C~ + + C . We see that the perturbative result is more than 3 times 
as large as the covariant one. Notice also that the perturbative results would reproduce very 
well the experimental amplitude (|8|), though one should recall that Coulomb effects would 
suppress the cross section at threshold by about 40 %. 

We do not consider here the model-dependent pion-production mechanisms, e.g., the 
pion-rescattering, and the p—u—ir diagram introduced in Ref. ||. They can easily contribute 
by the same amount as the model independent nucleon Born diagrams considered so far. 
However, care should be taken to include them in a way consistent with the dynamics of the 
AW-interaction. 

In conclusion, the sensitive cross section of 7r° production seems to be an ideal place 
to look for effects of relativistic dynamics. Our exploratory studies reported in this letter 
indeed lead to an assessment of the importance of "Z-diagrams" different from the tradi- 
tional nonrelativistic ones. Our main conclusion is that such diagrams, when included in 



a non-perturbative way, contribute differently than their non-relativistic limits taken per- 
turbatively. A similar effect was reported in Ref. [|D]] for the elastic electron scattering on 
deuteron, but there it becomes noticeable at a momentum transfer ~ 1 GeV. Our present 
calculations in relativistic impulse approximation underestimate the data. In future calcu- 
lations we plan to go beyond the threshold approximations and include pion rescattering 
contributions. 
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FIG. 1. The relativistic impulse approximation is denned by diagrams (a), (b), and (c). The 
shaded areas are NN scattering amplitudes. Diagram (d) is an example of the model dependent 
diagrams not included here. The thick lines represent off-shell nucleons, they can propagate with 
positive and negative energy. 
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FIG. 2. Total cross section vs lab energy of incoming proton. The Indiana data (]1| are shown 
with the best fits using unmodified p{\) and modified p(£) phase space densities, respectively. The 
curves labeled v = 1.6 and v = 0.0 are the full results of the covariant calculations. For the model 
with v = 0.0 we also show the cross section calculated with A4 + only, labeled " v = 0.0 +", to be 
compared with the perturbative IA "i/ = 0.0, pert IA" . The full perturbative prediction "i/ = 0.0, 
pert" virtually coincides with the p(£) fit. 
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FIG. 3. Ratio of half off-shell and on-shell amplitudes of the model with v 
two-nucleon lab energies (in MeV) as a function of the off-shell momentum p. 
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TABLES 

TABLE I. Amplitudes of diagrams (a), (b), (c) for several AW-interaction models. For dia- 
grams (a) and (c) the amplitude is also separated into contributions due to positive and negative 
energy nucleons in the intermediate state(s). To diagram (b) at threshold only intermediate neg- 
ative energy state contributes. M + = A + + C++, M~ = A~ + B~ + C+~ + C~+ + C . The 
correspondence between the real numbers A, B, and C and the complex amplitudes of diagrams 
(a), (b), and (c), is explained in the text. The label "no N ffs" refers to the total covariant results 
with the off-shell nucleon form factors switched off. Perturbative results are denoted by "Pert.". 
All amplitudes are given in units of 10 _7 MeV . 



Diagram 


v = 


v= 1.0 


u= 1.6 


i/ = 2.0 


A+ 


1.11 


-0.70 


-1.08 


-1.28 


A~ 


-1.44 


-0.87 


-0.04 


0.68 


A 


-0.33 


-1.58 


-1.12 


-0.60 


B = B~ 


-0.47 


-0.55 


-0.55 


-0.54 



C++ -6.28 -6.06 -5.73 -5.62 

C+- -5.42 -6.50 -6.41 -6.28 

C~+ 0.56 0.42 0.19 0.02 

C— 0.06 0.06 0.04 0.01 

C -11.1 -12.1 -11.9 -11.9 



-6.28 


-6.06 


-5.73 


-5.42 


-6.50 


-6.41 


0.56 


0.42 


0.19 


0.06 


0.06 


0.04 


-11.1 


-12.1 


-11.9 


-5.17 


-6.77 


-6.81 


-6.71 


-7.44 


-6.77 


11.88 


-14.21 


-13.57 


12.64 


-14.70 


-13.91 



M + -5.17 -6.77 -6.81 -6.90 

Mr -6.71 -7.44 -6.77 -6.11 

Total -11.88 -14.21 -13.57 -13.00 

no N ffs -12.64 -14.70 -13.91 -13.30 



Pert. I A 


-7.31 


-8.72 


-8.50 


-8.48 


Pert. Z-diag 


-23.3 


-24.6 


-24.8 


-23.7 


Pert, total 


-30.6 


-33.3 


-33.3 


-32.2 
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